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Abstract 

A new parameter, here termed openness, expressing the degree 
of dominance or enclosure of a location on an irregular 
surface, is developed to visualize topographic character. 
Openness is an angular measure of the relation between 
surface relief and horizontal distance. For angles less than 90°, 
it is equivalent to the internal angle of a cone, its apex at a 
DEM location, constrained by neighboring elevations within a 
specified radial distance. Openness incorporates the terrain 
line-of-sight, or viewshed, concept and is calculated from 
multiple zenith and nadir angles — here along eight azimuths. 
Openness has two viewer perspectives. Positive values, 
expressing openness above the surface, are high for convex 
forms, whereas negative values describe this attribute below 
the surface and are high for concave forms. Openness values 
are mapped by gray-scale tones. The emphasis of terrain 
convexity and concavity in openness maps facilitates the 
interpretation of landforms on the Earth's surface and its 
seafloor, and on the planets, as well as features on any 
irregular surface — such as those generated by industrial 
procedures. 

Introduction 

Relief maps and the parametric characterization of topography 
are essential to the interpretation of the land surface. Hill- 
shaded maps initially were drafted by landscape artists, but 
map precision and reproducibility were limited by the manual 
technique and its requisite cartographic skill (Raisz, 1931; 
Imhof, 1965; Alpha and Winter, 1971). The advent of the elec- 
tronic computer and the digital elevation model (DEM) over- 
came these restrictions, and machine-made shaded-relief 
maps now are widely used to display topography (Yoeli, 1967; 
Horn and Brooks, 1989; Thelin and Pike, 1991; Vigil et al., 
2000). The technique is not without drawbacks of its own. 
Because a directional light source is required, for example, 
ridges and valleys intersecting that source are shown clearly 
while features parallel to the source can be difficult to identify. 
Attempts to mitigate this problem by multiple light sources 
have been only partially successful (Mark, 1992; Moore and 
Mark, 1992; Riehle etal., 1997). 

Other types of neighborhood operations centered on DEM 
grid points yield digital maps that display or extract topo- 
graphic features (Tobler, 1969; Peucker and Douglas, 1975; 
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Weibel and Heller, 1991; Blaszczynski 1997; Guth, 2001). Ter- 
rain slope, curvature in both the XYand Z domains, relative 
relief, and slope azimuth are among the measures most com- 
monly mapped from DEMs (Evans, 1972; Pike, 1988). Algo- 
rithms that implement the much-used DEM-to-watershed 
transformation, now an essential tool for hydrologic and geo- 
morphic modeling, incorporate several such measures (Jen- 
son and Domingue, 1988; ESRI, 1992). Like automated relief- 
shading, however, all of these parameters are sensitive to the 
vagaries inherent in gridded height data (Carter, 1992). In par- 
ticular, slope maps computed from DEMs commonly have an 
unsatisfactory wormlike appearance that reflects the method 
of contour gridding and obscures the true configuration of the 
land surface. 

In this paper, we introduce a new concept of surface repre- 
sentation, provisionally termed openness for want of a better 
name. It is a novel method for digital terrain modeling — 
actually an image-processing technique— by which values of 
surface openness are calculated from a DEM, displayed in map 
form, and used to visualize landscapes. The resulting maps of 
openness superficially resemble digital images of shaded 
relief or slope angle, but emphasize dominant surface concav- 
ities and convexities. Values of openness require no light 
source, thus removing one limitation of relief shading, and 
are less affected by the DEM noise that afflicts most other 
parameters. After describing the technique, we present sev- 
eral maps of surface openness and discuss how they may be 
applied to the analysis and interpretation of topography. 

Zenith and Nadir Angles 

Before defining topographic openness, we must establish some 
angular relations between the point locations along individual 
DEM-derived profiles, and do this for viewer perspectives both 
above and below the surface. The first quantity to be extracted 
from these relations is elevation angle 0. 

Figure 1 shows a sample of terrain heights in a DEM arrayed 
north-south and east-west at a constant spacing Af, projected in 
Universal Transverse Mercator (UTM) coordinates, and cen- 
tered on a point of interest (double circle). Each point in the 
DEM is described by (i, ;, H) where i and / are column and row 
numbers and H is elevation, Figure 2 defines the geometric 
relation in profile between any two points A(i A , / A , H A ) and B(i B , 
/ B , H b ) in the DEM, here showing a case where H A < H B . Hori- 
zontal distance P between A and B is given by 
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Figure 1. Selected elevations arrayed in a DEM. Dou- 
ble circle is origin point for the calculation of open- 
ness; solid points are heights outward from that point 
along each of the eight azimuths; azimuth angle D 
is measured clockwise from north. L is radial limit of 
the calculation. 
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Figure 2. Elevation angle 6 between two points on a dem. 
Terrain height A(f A , / A , H A ) is at the origin of the calcula- 
tion; height B(/ B , J B , H b ) lies at some horizontal distance 
P on the dem. 9 is assumed to be positive for H A < H B , 
and negative for H B < H A . 



P = MjV A -i s r + lj A -h) 2 (U 
and the elevation angle 0 of line AB is 

Angle 6 is positive when H A < H B and negative when H B < 
H A . A general treatment of P and 6 would incorporate Earth's 
curvature, but this can be neglected because the calculation 
requires only small patches of terrain and thus short horizontal 
distances. Angle 0 enables us to identify aD - L set. 

Definition l: A D - L set for a DEM grid point is the set of all eleva- 
tion angles (0) between that point and each of the grid points 
(filled circles in Figure 1) located on a profile along an azimuth 
D and within radial distance L. The largest 6 in the D - L set, 
shown in Figure 3, is required to calculate elevation angles. By 
expressing the D - L set as D S L , we introduce two quantities, 
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Figure 3. Surface openness defined in terms of zenith and 
nadir angles. Zenith angle D <f> L or (90 - a f} L ), and nadir angle 
D iff L or (90 + „$.), calculated along one of eight azimuths 
D within radial limit L in Figure 1; dots are height along 
terrain profile. Positive openness is the mean value of D <fr. 
along the eight sampling directions; negative openness is 
the corresponding mean of eight values of D i[r L . Elevation 
angles ©A and d$. can be positive or negative, depending 
on the character of the topography around the central 
point A. 



n /3 L : the maximum elevation angle in the D S t 

and 

0 S L : the minimum elevation angle in the D S L . 

These angles, the minimum angles of elevation or depres- 
sion, respectively, for which the line-of-sight is unobstructed 
out to a specified range (here L), are the mask angles of the mili- 
tary tactician (Wood, 1963; Anon., 1990). These two angles, 
which also are involved in the calculation of terrain viewsheds 
(Lee, 1994, Fisher, 1996), enable us to define the next quantity. 

Definition 2: The zenith angle at a DEM grid point along azimuth 
D within radial distance L is 

D<k = 90 - D p L (3) 

and the nadir angle 

D <fr L = 90 + D S L (4) 

As illustrated in Figure 3, the zenith angle is the maxi- 
mum vertical angle subtended by a selected grid point and any 
of the points viewed from above the surface along a chosen azi- 
muth D up to distance L. The corresponding nadir angle D V h 
for that point is the maximum angle subtended by the chosen 
point and any other point viewed from beiowthe surface along 
azimuth D within L. Like terrain mask-angles , and D W L thus 
depend both upon the configuration of the surface sur- 
rounding the point of interest and upon the specified range L 
within which points are considered. Moreover, because D S L 
computed within increasingly large distances L includes more 
grid points, but not necessarily a commensurate increase in 
terrain relief, both ufa and D 1f £ tend to diminish with increas- 
ing L. 

Openness 

We obtain surface openness by calculating either zenith or 
nadir angles for all eight compass directions from the central 
point A (double circled in Figure 1) and taking the mean. 
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Definition 3: Positive openness &l at a location on a si 
within the distance of L on the DEM is 



and negative openness ¥ L within a distance L is 

= Wl. + 45& + ... + 315&)/8. [6) 

Positive openness at a DEM grid point A is the average of 
zenith angles subtended by A and the optimal point (that yield- 
ing the highest angle] viewed above the surface along each of 
the eight azimuths D 0 . 3 i5 — U P t0 ^ e specified distance L (Fig- 
ure 3). Similarly, negative openness averages the eight nadir 
angles subtended by point A and the optimal point viewed 
below Hie surface along azimuths Do-315 within L. To forestall 
possible confusion, we note that values of 4^ and ¥ L always 
take the positive sign. Openness is designated "positive" and 
"negative" to avoid the awkward terminology cited in the next 
paragraph. We adopt this "positive'V'negative" convention in 
the same sense as has been used to express terrain-slope curva- 
ture (Pike, 1988), which also takes only positive values (posi- 
tive curvature is convex-upward; negative curvature is 
concave-upward) . 

Positive, or "above-ground," openness <P L at a grid point is 
constrained by surrounding topography and thus is an angular 
estimate of the 360° horizontal extent of area around the point. 
We originally conceptualized the measure as "openness of the 
terrain to the sky." Its negative counterpart, "below-ground" 
openness W L — less helpfully verbalized as openness per se — is 
the corresponding angular measure of 360° spatial extent arro- 
und the central grid point. Figure 4, which shows these two 
measures schematically for angles < 90°, represents openness 
as the largest right cone of specified base diameter 2L that can 
be fit to a point on a DEM. Comparable representations of <Pl and 
V L > 90° are impossible-to-illustrate constructs in topology 
(e.g., Barr, 1964) that are better rendered in two dimensions, as 
in Figures 5A and 5D. Figure 4 is a simplification of the true 
figure of merit, because the base of the conical figure actually fit 
to the terrain rarely will be a perfect circle of diameter 2L. As 
shown in Figure 6, the base will have an irregular, commonly 
elongate and eccentric, form that necessarily varies with the 
elevation angles determined along the eight azimuths. 

Like its constituent individual zenith (d<&) and nadir Cd^J 
angles, the corresponding mean values of openness and ¥ L 
also are constrained to decrease with increasing range L. 
Accordingly, the choice of L will emphasize topography at dif- 
ferent length scales. Larger values of L will highlight larger fea- 
tures, and smaller L smaller forms. 

Discussion 

Openness expresses a topographic attribute that possibly is 
captured by no other measure or combination of them, 



although this remains to be tested against parameters that 
describe terrain fabric (Pike et ah, 1989; Guth, 2001). The char- 
acteristic appears to be the degree of geometric "dominance" of 
one location on a surface by another. This parameter expresses 
two configurations commonly encountered in the landscape: 
command of an expanse of terrain by one or several elevated 
relief features (positive openness and the degree of "enclo- 
sure" of a lower location by elevated surroundings (negative 
openness Sy . As such, openness is sensitive to the local relief ; 
distance conditions that express much of an area's topographic 
distinctiveness (Pike etal., 1989). 

The range of contrasts is illustrated by profiles of a hilltop 
and a closed depression along the azimuth D w -D 270 (Figure 5). 
If the center point of an openness calculation is on a summit 
rising above its surroundings and flanked by steep slopes, the 
resulting values of 90 <£t and 2 7o<k. will exceed 90° (Figure 5 A), as 
will positive openness , which is the average of eight values 
of d^l. S ummi ts and ridge lines thus will appear on a plot of <Pl 
as, respectively, compact and linear groupings of high-scoring 
cells. Commensurately, a negative-openness calculation at the 
same location (Figure 5B) would yield low values of x i)t L and 
270^1. < 90°, and summits and ridge lines on a plot of 9% will be 
groups of low-scoring cells. Topographic features that result in 
neither high nor low values of <Pj, and V L yield groupings of 
intermediate-scoring cells. Taking the inverse example, if the 
center point of an openness calculation lies within a depression 
surrounded by steep slopes (Figure 5C), the resulting value of 
positive openness &l will be less than 90 degrees. Topographic 
lows, therefore, will appear on maps of & L as groups of low- 
scoring cells. However, that same depression will generate a 
high value of negative openness 9i > 90° (Figure 5D). Topo- 
graphic depressions and valley lines thus will show up on a (L 
plot as clusters and strings of high-scoring cells. 

The closed-depression example, which is the approximate 
inverse of the summit example, raises an intriguing point. Fig- 
ure 5 would at first appear to show that negative openness is 
redundant: that simply subtracting <t> L from 180° would be 
equivalent to but this is not so. The two quantities differ 
substantially, due to irregularity of the ground surface sur- 
rounding point A. It is for this reason that negative openness V L 
emphasizes concave features (Figure 5D) more crisply in gray- 
tone images (Figures 7 through 13 below) than does 180° - S L 
(Figure 5C). Otherwise, 180° - would suffice to capture 
depressions and valleys and there would be no need for Tfi,. 

The correspondence of openness to terrain configuration is 
illustrated in Figure 6 by octagonal plots that compare zenith 
and nadir angles, along the eight azimuths, for nine test sur- 
faces. Among these simple models are all the important ele- 
ments of topography — plain, slope, summit, depression, ridge, 
valley, and saddle (Maxwell, 1870). Depending on shape of the 
ground, angular values of D <f> L and D ifi L increase with distance 
outward from the center of each figure, up to range L. Openness 
is equivalent to the area of the octagon projected at L. Because 





Figure 4. Positive (left) and negative (right) openness shown schematically for 
values < 90°. Heavy irregular line is terrain surface; L is radial limit of calculation 
for chosen point (large dot) on a dem. 
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Rgure 5. High and low values of positive # t and negative ¥ L surface open- 
ness, illustrated. Profiles of a hilltop and closed depression along azimuth 
D90-D270 (E-W). Values of 30 <f> L and 2 io<k for point Aj. each > 90° in (A), but 
values of 90 i/r L and 27 o& < 90° In (B). Because <ti and V L are means of 
D 4> L and D tj/ L along eight azimuths, 4> L at A x > 90° and 9i < 90°. Similarly, 
90 <k and 27o<k for point A 2 each < 90° in (C), but those of 90 & and 270& 
> 90° in (D), thus <P L at A 2 < 90° and > 90°. None of these angles is 
larger 180°. 



features in natural terrain are complex composites of these 
basic elements, the set of eight angles that comprise openness 
will vary considerably with location of the point of origin and 
character of the surrounding topography. A grid point in highly 
irregular terrain thus will generate a variety of constituent 
angles and a correspondingly distorted octagonal pattern. The 
angular plots in Figure 6 constitute a kind of geometric signa- 
ture of topographic form; openness also could be part of a more 
complex, multivariate signature (Pike, 1988). 

Examples and Applications 

While we defer the application of corresponding values of 
openness as an analytical parameter — for example, as a con- 
stituent of a geometric signature — to a later experiment, in this 
paper we exemplify the visuahzation counterpart of terrain 
modeling by several maps of surface openness. These few selec- 
tions represent a variety of terrains, map scales, L-values, and 
DEM sources and resolutions. Not all possible combinations of 
these applications are shown, hi order to highlight both terrain 
convexities and concavities, for example, values of positive 
openness and negative openness must be compiled as separate 
maps of the same area. Also, separate openness maps of the 
same area can be computed for different values of L, to empha- 
size fine- or coarse-scale features . 

As in mapping elevation, slope, and shaded relief calcu- 
lated from DEMs, tones of gray are assigned to values of open- 
ness within a range that produces an image with optimal detail 
and contrast. The convention adopted for our images is that 
customarily followed for elevation and slope: high values of 



3 and low values by dark 
tones. Accordingly, images of positive openness <P L , which are 
designed to highlight topographic convexities, show ridges as 
bright lines and rivers and other elongate concavities as dark 
lines, whereas maps of negative openness V L emphasize drain- 
ages (light tones) at the expense of mountains and other convex- 
overall features (dark tones). Intermediate values of both 
parameters yield gray tones. Although we experimented with 
reversing this tonal convention, we chose the standard scheme 
on aesthetic and intuitive grounds — the resulting images sim- 
ply look more pleasing to the eye and seem easier to interpret. 
Mapping openness values carries some perceptual conse- 
quences. While openness images superficially resemble raster- 
based maps of elevation, shaded relief, or slope, their appear- 
ance can be deceiving at first. Individuals familiar with the 
more conventional types of digital maps may require some 
acclimation to openness images before viewing them becomes 
intuitive. 

Maps of positive openness in Figures 7 and 8, computed at 
different values oft, illustrate various features of stratovol- 
canoes on the main Japanese island of Honshu (Yoshikawa et 
al., 1981, p. 81-87). Landforms in the vicinity of Mt. Fuji in 
south-central Honshu (Figure 7) that are easily recognized from 
openness textures include aligned low hummocky hills 
formed in flank eruptions prior to 1707 (light-toned area, cen- 
ter). Steep ridges in much of the image appear as bright thin 
lines; old lava flows on Fuji and the more gently sloping ridges 
on the eroded Mt. Ashitaka volcano and Hakone caldera to the 
Southeast are less distinct. The area around Mt. Chokai in 
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Figure 6. Variation in surface openness. Differing values of the eight zenith angles D <pi and nadir angles D i]/ L , at 
a dem grid point (black dot), for nine simple terrains. Positive and negative openness <P L and tfi. are proportional 
to the area enclosed by the octagonal plots of D <k and D A, respectively. Evenly-spaced numbers 1 through 7 
along azimuths are for comparison purposes and do not correspond to actual angles; range L for all calculations 
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Figure 7. Positive openness (L = 5 km) for the Mt. Fuji area, southern 
Honshu <N35:06:24 - N35:31:42, E138:18:07 - E139:09:16) com- 
puted from a 50-m dem of the Geographical Survey Institute of 
Japan (gsij). 




Figure 8. Positive openness (1 = 1 km) for Mt. Chokai, 
northern Honshu (N38:43:41 - N38:59:32, E139:44:49 - 
E140:14:56) computed from a 10-m dem from the Hokkaido 
Map Company (hmc). 



northwest Honshu (Figure 8), rendered at one-fifth, the value of 
L in Figure 7, is more detailed and also differs in morphology 
from the area around Mt. Fuji. The recent growth process of the 
volcano, clearly evident in the lobate form of its many fresh lava 
flows, is well illustrated by positive openness. 

Figures 9 and 10 contrast two non-volcanic terrains in 
Japan in maps of negative openness that emphasize topo- 
graphic concavities at very different values of L. The Tohka- 
machi area (Figure 9), located in the upper reaches of the Shi- 
nano River in north-central Honshu, exemplifies variation in 
the patterns of valleys as well as in the layers of sediment in the 
main valley. Figure 10 shows the Tohno area, east of Mt. Chokai. 
The distinctly reticular structure of this terrain, at one-tenth 
the L used to compute Figure 9, reflects small stream channels 
that developed along planes of weakness caused by jointing in 
the underlying granite bedrock Because the area is densely 
vegetated, the macroscopic structure of this joint pattern had 
remained unknown to geologists and is shown clearly here for 
the first time (Kanisawa and Yokoyama, 1999). More variety of 
topographic features can be observed in the openness maps of 
Japan (Yokoyama, 2001), which were generated from a 50-m 
interval DEM with L = 5 km and printed at a scale of 1:50,000. 



Figures 11 and 12, showing positive and negative openness 
at the same L-value for a small area 50 km south of San Fran- 
cisco in California's coastal Santa Cruz Mountains, are 
excerpted from two large openness maps of the entire ten- 
county San Francisco Bay region (not shown). Most of this 
rough terrain lies within the La Honda 7.5-minute quadrangle, 
the first topographic map published by the U.S. Geological Sur- 
vey in which a contour-to-grid technique replaced photogram- 
metric scanning to create a higher quality 30-m DEM (Pike et al., 
1987). The crisp valleys and ridge lines in Figures 11 and 12, 
respectively, reflect the region's strong contrasts in rock type 
(Pike ef al. , 1988) as well as the imprint of recent tectonic uplift 
and consequent erosion. The steep terrain is subject to chronic 
debris-flow (landslide) activity during the rainy California 
winter and thus is thinly populated. 

The lofty topography of Tibet and the high-relief Himalaya 
Range are attracting fresh scientific interest (Fielding ef al., 
1994; Albright ef al., 1998; Zeitler ef al., 2001). Figure 13, aneg- 
ative openness map of this large region, emphasizes the 
hydrology of the southern Tibetan Plateau and its relation to the 
Himalayas. The distinctive north-south valleys were carved 
through the rising mountains by rivers draining the plateau. 
The contrasting rectangular areas reflect differences in the 
source and quality of data in the GTOP03D DEM. This contrast is 
an important feature of the openness parameter. The large map 
of negative openness (not shown) containing the California 
sample in Figures 11 and 12 also revealed an imperfect match 
between two source DEMs — suggesting that the openness is an 
excellent tool for performing quality control on regional DEMs 
mosaicked from several smaller data sets. 

Application of openness need not be limited to the por- 
trayal and interpretation of Earth's landscapes. Like shaded- 
relief mapping, openness mapping is an image-processing tech- 
nique capable of representing surface features anywhere aDEM 
exists or can be prepared — Earth's seafloor (Smith and Sand- 
well, 1997) and surfaces of the solid planets and satellites 
(Smith et al., 1999). Finally, the technique is not restricted to 
natural landscapes. Digital terrain modeling has an important 
micro- and nano-scale counterpart in surface metrology, the 
numerical characterization of industrial surfaces — such as 
magnetic tape and disk surfaces and automobile-engine cylin- 
der walls (Pike, 2000; Pike, 2001). Metrology recently was revo- 
lutionized by its ability to create and manipulate square-grid 
height matrices that are fine-scale analogs of DEMs. Because the 
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Figure 9. Negative openness (L = 5 km) for the Tohka-machi area, 
central Honshu (N36:54:07 - N37:19:20, E138:21:58-E139:09:05) 
computed from the 50-m gsij dem. 




Figure 10. Negative openness (t = 500 m) for Tohno area, 
central Honshu (N39:ll:00 - N39:36:41, E141:09:58 - 
E142:01:06) computed from a 10-m dem from hmc. 




Figure 11. Positive openness (L = 500 m) for an area south 
of San Francisco, California (N37:10:03 - N37:17:43, 
W122:06:01 - W122:22:l) computed from a 30-m dem of 
the U.S. Geological Survey (usgs). 




Figure 12. Negative openness (JL = 500 m) for the area in 
Rgure 11. 



surface topography of manufactured components is routinely 
modeled by the same techniques of parameterization and visu- 
alization as used in the Earth sciences, we suggest that the con- 
cept of surface openness offers a new tool for metrologic 
applications as well. 

Conclusion 

A new image-processing technique generates an angular meas- 
ure of surface form, here termed openness, that visualizes the 
topographic dominance or enclosure of any location on an 
irregular surface represented by a DEM. The measure incorpo- 
rates the terrain Iine-of-sight (viewshed) principle and is calcu- 
lated from zenith and nadir angles along eight DEM azimuths. 
Openness is expressed in two modes. Positive openness 
emphasizes convex features of topography and negative open- 
ness concave features. Values of openness, which are shown as 
gray-scale maps resembling shaded-relief images, are indepen- 
dent of the position of any light source and less sensitive to 
noisy DEMs than are other parameters. Different topographic 
features can be emphasized by tailoring a radial sampling dis- 
tance to the scale of the selected area. Openness maps not only 
create new possibilities for geomorphological and geographical 
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Figure 13. Negative openness (L = 5 km) of the Tibetan Plateau and 
Himalaya (N25:06:30 - N38:14:00, E75:30:30 - E97:49:00) com- 
puted from a 1-km GTOP030 dem available from the usgs. 



interpretation, but also are suited to visualization of the Earth's 
seafloor, planetary landf orms, and features on any irregular sur- 
face, including those micro- and nano-landscapes created by 
manufacturing processes. 
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Representation of topographical features by opennesses 
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Abstract Overground-openness and underground-openness are defined for a grid 
5 point on digital elevation models. The former is a characteristic quantity to describe sky extent 
over the point within a distance L and takes large values for convex landforms. The latter is to 
describe underground extent and takes large values for concave landforms. These concepts are 
applied to 5Qm-mesh digital elevation model of Mt Iwate to make openness maps of L = 5 km 
and L = 0.5 km The maps were found useful to represent topographical features such as the 
1 0 lines of ridge and valley. 
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1. Introduction 

1 5 Diffusion of digital elevation model has allowed effective topographic representation 

and analysis by computer. As a method of displaying a topographic map with an appearance of 
solidity, there has been known from old times a technique of relief map by a hill shading method 
This is a method of indicating slopes facing a tight source, brightly, and indicating slopes with 
the light source in the background, darkly. Although manually produced relief maps were 

20 subject to limitations in accuracy and reproducibility, now with mathematically systematized 
procedures for preparation of a relief map from a digital elevation model, it is allowed to prepare 
various ones as necessary for the use 1 2) . In the relief map, changes in terrain crossing an arrival 
direction of light can be clear represented, but changes in terrain parallel to the arrival direction 
of light become unclear. On the other hand, there have been methods proposed for extraction of 

25 a variety of topographic features (ridgeline, valley line, gradient, valley density, water catchment 
area, etc.) from a digital elevation model However, most of their algorithms have been based 
on a local calculation by elevation values in a vicinity of a site, susceptible to noises, and unable 
to grasp comprehensive topographic features, as shortcomings. 

This paper introduces a concept of "openness", to propose a new technique for 

30 representation of topographic features based on a digital elevation model. Viscerally, the 



openness is a quantification of a degree by which the site is projected on the earth in comparison 
with peripheral, or a degree by which the site is cut belowground. Mapping the openness 
permits extraction of information suitable to a topographic scale depending on a specified 
calculation distance, allowing for a representation free of orientation or local noise. Actually, 
5 openness maps of areas around Mt Iwate were prepared on the basis of a "numerical map 
50m-mesh (elevation) " issued by the Geographical Survey Institute of Japan, revealing that the 
present technique is particularly excellent in extraction of ridgelines and valley lines, allowing 
wealth of topographic and geologic information to be made out 

There will be described in Chapter 2 the "numerical map 50m-mesh (elevation) ", and 
10 in Chapter 3, the openness will be defined. In Chapter 4, prepared openness maps will be 
discussed. Chapter 5 provides a conclusion. 

2. Preparation of a digital elevation model by UTM coordinate system 

The "numerical map 50m-mesh (elevation) " is a two-dimensional array of elevations, 

15 read at intervals of 10 cm, at centers of meshes (mesh intervals being 2.25 seconds in the 
latitudinal direction, and 1.50 seconds in the longitudinal direction) equally dividing the height 
and the width of a 1/25,000 terrain map. In the "numerical map 50m-mesh (elevation) ", the 
distance between sample points (as positions where elevations are read) varies depending on the 
longitude and latitude (for example, on a northern latitude line of 40 degrees, relative to the 

20 north-south direction where it is 46.26 m in, it becomes 53.37 m in the east-west direction), in 
addition to that arrayed directions of sample points are disaccord with azimuths, except for the 
north-south, such that the greater the distances the greater the differences. 

For a convenient topological analysis over a wide area, the digital elevation model 
should have sample points at a constant spacing, and data array directions in accord with 

25 azimuths. Hence, based on the "numerical map 50m-mesh (elevation) ", we prepared a set of 
digital elevation data by an UTM projection (in application of Zones 52 to 55) having a sample 
point spacing of 50 m This will be referred to as "UTM-DEM'. For conversion from the 
"numerical map 50m-mesh (elevation) " to the UTM-DEM, there was employed a conformal 
first-order interpolation 7 K 

30 Fig. 1 illustrates an arrangement of sample points in the UTM-DEM As the azimuth is 



defined in terms of a clockwise angle D relative to the north to be 0 degrees, there are eight 
azimuths (north, northeast, east, south east, south, southwest, west, north west) appearing in 
sequence with respect to each 45 degrees from 0 degrees. Hereinafter, a single sample point will 
be described by 
5 (U,H) 

where i and j are numbers of a row and a column indicating a position of the sample point, and 
His the elevation. 

10 ♦ 



15 




D=iao 



20 Fig. 1: Sample points arrayed about a focused sample point (with a double circle 

mark) in the UTM-DEM, in which black circle marks are sample points along the eight 
azimuths, and white circle marks are sample points else. 

3 . Overground openness and underground openness 
25 All the following discussions presuppose the UTM-DEM. 

3 . 1 Overground angle and underground angle 

Suppose two sample points A (ia, ]a, Ha) and B (i B , j B , H B ). For the sampling at 
intervals of 50 m, the A and B have a distance in between, such that 
30 P=50^l(i A -i B ) 2 +U A -JB) 2 




Fig. 2 illustrates a relationship between sample points A and B, with an elevation of 0 
m as a reference. The sample point B has an elevation angle 6 with respect to the sample 
point A, to be given such that 

0 = imx- l {(H B -H A )/P} 
5 where 9 has a sign to be: (D positive for H A < Hfe; or (D negative for H A >H B . 



Plane of a 0-m elevation 



Fig. 2: Elevation angle d of a sample point A with respect to a sample point B 

With respect to a focused sample point, a set of sample points along an azimuth D 
wilhin a range of a distance L is described by d Sl, which will be referred to as "D-L set of the 
1 5 focused sample poinf \ 

Now (referring to Fig. 3), letting 

D /3 L : the maximum among elevation angles of set elements of a d Sl of a focused 
sample point, and 

D 5 L : the minimum among elevation angles of set elements of a D S L of a focused 
20 sample point, 

a definition is made, as follows: 




(over the ground) 



Fig. 3: Overground and underground angles of a D-L set of a focused sample point A. 



Definition I: With respect to a D-L set of a focused sample point, an overground angle 
and an underground angle mean 
d <£l = 90-d£l, 

5 and 

d ¥ l = 90 + d 6l, 
respectively. 

D <j> L means the maximum of a zenithal angle by which a sky along an azimuth D is 
visible within a range of distance L from a focused sample point What 1hey call skyline angle 

10 corresponds to the zenithal angle, as the L is set infinite. And,^!. means the maximum of a 
nadir angle by which a nadir along an azimuth D is visible within a range of distance L from a 
focused sample point For an increased L, the number of sample points belonging to the d Sl is 
increased, so the D B l has a non-decreasing characteristic with respect to L, while the D 8 l has a 
non-increasing characteristic, vice versa. Accordingly, both d 4> l and d ^ l have a 

1 5 non-increasing characteristic with respect to L. 

In the surveying, the elevation angle is a concept denned to a horizontal plane as a 
reference passing a focused sample point and does not strictly accord with 9 . Further, for strict 
discussions of the overground and underground angles, also a curvature of the earth should be 
considered, so the definition I is not always an exact description. The definition I is a concept 

20 defined absolutely assuming the use of a UTM-DEM for topographic analysis. 

3 .2 Overground openness and underground openness 

The overground angle as well as the underground angle is a concept with respect to a 
specified azimuth D, which is now extended to introduce a definition, as follows. 
25 Definition II: For a distance L with respect to a focused sample point, an overground 

openness and an underground openness mean 

$ L = (o 4> L + 45 <2> L + 90 4> L + 135 4> L + 180 $ L + 225 L + 270 4> L + 3 15 <f> L ) / 8 

and 

1 FL=(o 1 FL + 45^L+90^L + 135^L + 180^L+225^L + 270^L+315^L)/8, 

30 respectively. 



The overground openness represents an extent of sky visible within a range of distance 
L from a focused sample point, and the underground openness represents an extent of nadir 
within a range of distance L„ when looking under the ground, standing on the head (refer to Fig. 
5 4 [overleaf]). 

(Underground openness) 
Fig. 4: Opennesses at a focused sample point (marked with a black circle) 

1 5 The opnennesses depend on the distance L and a peripheral terrain. Listed in Table 1 

are sets of azimuth-wise octagonal plots of overground angles and underground angles 
representing overground opnnesses and underground opennesses of nine fundamental 
topographic forms. Generally, the higher a site is projected from around the greater the 
overground openness becomes, so the value is large on a peak or ridge, and small in a dent or 

20 valley floor. To the contrary, the lower a site is cut underground the greater the underground 
openness becomes, so the value is large in a dent or valley floor, and small on a peak or ridge. 
Really, there are various mixtures of fundamental topographic forms even within a range of 
distance L, giving various combinations of values to sets of octagonal plots of overground 
angles and underground angles. 

25 

Table 1 : Overground opnnesses and underground opennesses of focused sample points 
(marked with black circles) on nine fundamental topographic forms, represented by sets of 
octagonal plots on azimuth-wise relative scales calibrated with five graduations to opennesses of 
overground and underground angles of flafland (to be 90 degrees). 

30 
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(Overground openness) 



Fundamental topographic forms 



Underground 



Flatland 



(middle) 



Mountaintop 



(very large) , 

7^ 



(very large) 



Ridge 

extending 

North-South 



(large) 



VaUey 



North-Sourth 



(large) 



Col between 
E-W valley and 
N-S ridge 



X 



(middle) 



(middle) 



Uniform slope 
facing East 



(middle) 



Boundary between 
western flafland and 
eastern down-slope 



(middle large} 



(middle small) 



Boundary between 
western flatland and 
eastern up-slope 



(middle small) 



4. Openness maps prepared 

For an area about Mt Iwate, between latitudes 39 degrees 36 minutes and 39 degrees 
57 minutes north and between longitudes 140 degrees 44 minutes and 141 degrees 23 minutes 
east, openness maps of L = 5 km and L = 0.5 km were prepared. An operational processing by 
5 Laplacian also was performed for a comparison purpose. 



4. 1 Topographic features of target area 

Fig. 5 is a chorography of the target area. Mt Iwate (2,038 m) is substantially centered, 
with Shizukuisbi Basin at the southern foot At the eastern foot there is a flatland extending as 
10 upstream part of Kitakarni River. At east of the flatland, there are closing mountains in western 
part of Kitakarni Highlands about Mt Himegami (1,124 m). Further, there is Ohu' Range of 
Mountains at the west of Mt Iwate, with volcanic terrain formed along a district between 
Hachimantai (1,613 m) and Akita-komagatake (1,637 m), and in the southern west part, with 
non-volcanic terrain spread about Gobanmori (1,048 m). 




Fig. 5: Target area of prepared openness maps (between latitudes 39 degrees 36 
minutes to 39 degrees 57 minutes north, and between longitudes 140 degrees 44 minutes to 141 
degrees 23 minutes east) 

30 



4.2 5 km overground openness 

Fig. 6 is an overground openness map for L = 5 km The greater the openness the more 
white the indication, whereby tops and ridgelines extending from them appear as white, and 
valley lines, as black, while, among others, flatlands and slopes or such appear as intermediate 
5 gray tones. The more white a ridgeline the higher the latitudes in comparison with surroundings 
(within a range of calculation distance L), and the thinner it is the steeper the topographic form. 
Likewise, the darker a valley line the deeper the valley, and the thinner it is the narrower the 
valley floor. For the former, as a typical, there is a ridgeline extending from the top of ML Iwate, 
and for the latter, there is Kuzuneta Ravin extending along a southwest foot of Mt Iwate. The 

10 district at the west side of Mt Iwate is a typical shield volcanic (aspite) district forming a gentle 
slope, which is seen as a black-and-white patched texture. On the other hand, districts of 
Gobanmori along the southwest edge as well as Kitakami Highlands at the east form textures 
with featuring black and white lines, whereby it is found that they constitute an erosional 
topography with intersecting ridgelines and valley lines. It can be seen that, for Gobanmori, with 

1 5 great and rough black-to-white contrasts, the landform should have actively dissected steep and 
large undulations, and for Kitakami Highlands, with small and fine contrasts in texture, the 
landform should have small undulations. Flatlands in upstream part of Kitakami River and 
Shizukuisi Basin are flat, and those districts appear gray as a whole, while terraces as well as 
erosional landforms can be read when observed into details. 

20 
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(openness) 

3 0 Fig. 6: Overground openness map for L = 5 km 




4.3 5 km underground openness 

Fig. 7 is an underground openness map for L = 5 km. The greater the openness the 
more white the indication. The deeper a valley the more white it is, and the wider the floor the 
bigger the valley appear. Even in those parts extending as relatively white regions, more white 
5 lines are visible, which are valley floors with advancing new dissections. Overground 
opennesses in Fig. 6 and underground opennesses in Fig. 7 have correspondence relationships, 
though representing topologic features from different points of observation, respectively. From 
the viewpoint that they are both easy to make out by white parts to be focused, we can say that 
Fig. 6 is suitable to make out a topography mainly composed of ridgelines, and Fig. 7, suitable 

10 to make out a topography mainly composed of valley lines. It may well be said that important 
topographic and geographic information on undulation, dissection, swelling, and the like are 
reflected in black-and-white textures of Fig. 6 and Fig. 7, i.e., in parts of lineated texture, such as, 
for lines, on the density, whether long or short, thick or thin, the degree of flexion, or the contrast, 
and in parts of patched texture, such as on the contrast, extent, or relationship of connection. 

15 Fig. 6 and Fig. 7 exhibit comprehensive topographic features within a range of 5 km 

from a focused point, and are free from such defects as weak to noises of local calculations. 
Further, the openness is calculated as an average of overground angles or underground angles in 
eight azimuths, and substantially orientation-independent 




in 



4.4 0.5 km overground opennesss 

Fig. 8 is an overground openness map for L = 0.5 km Though being alike with Fig. 6 
as an overground openness, Fig. 8 is stronger in black-to-white contrast, with more appearing 
black and white lines. Since the openness has a non-increasing characteristic with respect to L as 
5 described, a topographic feature of a closer region to a focused sample point is to be extracted, 
as a smaller L is designated. Those ridgelines and valley lines appearing in Fig. 6 will appear in 
Fig. 8 as well but all ridgelines and valley lines in Fig. 8 will not appear in Fig. 6. In this sense, 
macroscopic topographic features are extracted in Fig. 6, and microscopic topographic features 
are extracted in Fig. 8, instead. 

10 
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Fig. 8: Overground openness map for L = 0.5 km 



4.5 0.5 km underground openness 
25 Fig. 9 is an underground openness map for L = 0.5 km In this case also, due to the 

difference in calculation distance L, relative to Fig. 7., the black-to-white contrast is stronger, 
with more appearing valley lines. In particular, among others, in foot parts of Mt Iwate as well 
as volcanic districts of Akita Komagatake and Hachimantai at the west side of Mt Iwate, minute 
valley lines can be seen, which are unobservable in Fig. 7. 

30 



1 1 



5 



10 




(openness) 



Fig. 9: Underground openness map for L = 0.5 km 

4.6 Topographic map by Laplacian 
15 The UTM-DEM of target area was processed by Laplacian operation, with a result 

shown in Fig. 10. The Laplacian is a known image processing measure for edge detection based 
on a spatial secondary differentiation * 
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(Laplacian) 

30 Fig. 10: Result of Laplacian operation processing 
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Generally, the Laplacian provides a large positive value to convex terrain, and a small 
negative value to concave terrain. Further, it gives near-zero values to landforms with small 
gradient variations, such as a flatland or monotone slope. In representation in the figure, the 
larger a Laplacian value is the more black the site becomes, and the smaller it is the more white 
5 the site becomes. Though actually appearing as black and white lines, the ridgelines and valley 
lines are matte, and have poor information readable, in comparison with Figs. 6 to 9. 



25 



5. Conclusion 

In this paper we have proposed, as new measures for representation of topographic 
features, an overground openness and an underground openness under a designated distance for 
calculation using a digital elevation model. Actually, for a peripheral area of Mt Iwate, we 
prepared overground openness maps and underground openness maps with different calculation 
distances. The overground openness maps and underground openness maps adequately reflected 
features of convex landforms and concave landforms, respectively. It also has turned out that 
macroscopic topographic features can be extracted by adoption of a large calculation distance, 
and microscopic topographic features can be extracted by adoption of a small calculation 
distance. In comparison with conventional relief maps or maps by differential operations for 
such subjects, the above subject maps contain rich topographic information, with possible 
expectation to open a new frontier for topographic and geographic decipher. Now, we are 
preparing overground openness maps and underground openness maps of extended target areas, 
with confirmed appearance of featuring images in figures, covering terrace, volcano, alluvial fan, 
hill side, faulting, and other landforms, as well. The opennesses have abilities to be effectively 
employed for topographic decipher and classification, as well, so we are repeating discussions 
with experts in associated fields. We want to report results of such researches another time. 
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